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We study the role of multiplicative colored noise for different values of the correlation time Tc 
in the dynamics of two competing species, described by generalized Lotka-Volterra equations. The 
^ ' multiplicative colored noise models the interaction between the species and the environment. The 

, interaction parameter between the species is a random process which obeys a stochastic differential 

' equation with a generalized bistable potential in the presence of a periodic driving term, which 

\l ' accounts for the environment temperature variation. The bistable potential is useful to describe the 

coexistence and exclusion dynamical regimes of the ecosystem. Noise-induced periodic oscillations of 
the species concentrations and stochastic resonance phenomenon appear due to the presence of the 
fH ^ multiplicative noise. We find that for low values of the correlation time Tc the response of the system 

■ coincides with that obtained with multiplicative white noise. For higher values of Tc the coherent 

D ' response of the system and the maximum of the signal-to-noise ratio, signature of the stochastic 

resonance phenomenon, are shifted towards higher values of the noise intensity. 
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I. INTRODUCTION 



Recently the constructive role played by the noise in population dynamics has been the subject of intense theo- 
retical investigations The noise, through its interaction with the nonlinearity of the ecosystems, can give rise to 
' ■ counterintuitive phenomena like stochastic resonance^ (SR), noise enhanced stability^ (NES), noise delayed extinction 

'S: ^etc...^ 

Q . In this work we focus on the role of realistic noise in the dynamics of two competing species, and specifically on 
^ ' the SR phenomenon in population dynamics in the presence of exponentially correlated noise. When the time scale 
of random fluctuations is larger than the characteristic time scale of the ecosystem the external noise cannot be 
^ ' considered white noise. This noise has no time scale because is J — correlated. In real ecosystems the external random 
^ ' perturbations, due to interaction with the environment, are correlated with a finite correlation time. A strongly 
[ — ■ correlated noise for example emerges as the result of a coarse graining over a hidden set of slow variables®. The 
00 dynamics of our ecosystem is described by generalized Lotka-Volterra equations with a random interaction parameter 
between two competing species in the presence of a multiplicative colored noise. In the absence of noise we have two 
stable states, corresponding to dynamical regimes of coexistence and exclusion. Real ecosystem however are subjected 
to noisy environment and to geological periodical forcing, like seasonal temperature. Because of these random and 
periodical forcing we describe the two dynamical regimes of coexistence and exclusion with a bistable potential for 
the interaction parameter, tilted periodically by an environment temperature driving term. 
' We find that the noise is responsible for the generation of quasi-periodic temporal oscillations and the enhancement 
^ of the response of the system to a driving force producing stochastic resonance^. By analyzing the effect of the colored 
I . noise we find that for low values of the correlation time Tc the response of the system coincides with that obtained 
''O ' with multiplicative white noise. For higher values of Tc the coherent response of the system and the maximum of the 
Ch signal-to-noisc ratio (SNR), which are signature of the SR phenomenon, are shifted towards higher values of the noise 
intensity. These results are in agreement with previous theoretical and experimental investigations of SR phenomenon 
in dynamical systems in the presence of colored noise®'^'^*^. However in previous studies the colored noise was additive, 
while here we have two different sources of noise and only one of them is colored. 



■ II. THE MODEL 

Time evolution of two competing species is obtained within the formalism of the Lotka-Volterra equations^ ^ in the 
presence of a multiplicative noise 

dx 

— ^ ^ii X {ai - X - l3i{t)y) + X Cx{t) (1) 
at 



dt 



— = ti2y{a2- V - P2{t)x) +yCy(t), (2) 
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where (i{t) {i = x, y) are colored noises given by the archetypal source for colored noise, i. e. exponentially correlated 
processes given by Ornstein-Uhlenbeck process^^ 

\i + -m {i = x,y) (3) 



dt 



Tn Tn 



and S,i{t) {i = x,y) are Gaussian white noises within the Ito scheme with zero mean and correlation function 
{Ci{t)^j{t')) = 2c7(5(f — t')Sij. The correlation function of the processes of Eq.(3) is 

(0(i)0(i')> = fe-l*-*'l/-<5,, (4) 

and gives 25{t — t')Sij in the limit 0. Wc can distinguish two different regimes for the correlation time of the 

noise: (a) dt < Tc < To, i. e. a correlation time greater then the integration step dt and less then the period of the 
deterministic driving force Tq; (b) Tc > Tq. 

The time series for the two populations arc obtained setting ai = a2 = a, l3i{t) = f32{t) = P{t), and £,x{t) = £.ij{t) = 
£,{t). For /? < 1 both species survives, that is a coexistence regime takes place, while for /3 > 1 an exclusion regime 
is established and one of the two species extinguishes after a certain time. Coexistence and exclusion of one of the 
two species correspond to stable states of the Lotka-Volterra's deterministic model^^. This aspect together with the 
assumption of a strong noisy interaction between species and environment suggests that the interaction parameter 
(3{t) can be described as a stochastic process which obeys the following stochastic differential equation 

dm _ du{p,t) ^ 

-IT- d^^^'^'^' 

where f/(/3, t) 

U{P, t) = h{(3 - (1 + p))V??' - 2/i(/3 - (1 + p))V^' - Picos{ivot), (6) 

is a bistable potential (see Fig.l), with the two stable states of coexistence and exclusion, modulated by the periodic 
term jcos^uJot), which takes into account for the environment temperature oscillations. Here h is the height of the 
potential barrier. 




FIG. 1. The bistable potential U{l3,t) of the interaction parameter P{t) at uiot = 7r/2. The potential U{(5,t) is centered on 
P = 0.99. The parameters of the potential are h = 6.25 • 10~^, r] = 0.05, p = —0.01. 



In Eq. (5) is a Gaussian white noise with the usual statistical properties: (Cpit)) ~ and {£,fj{t)£^fj{t')) = 

a5{t — t'). We choose as initial conditions for two competing species a coexistence regime (/3(0) < 1) on the left 
minimum of the potential U{f3,t) at io^t = it/2. 
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A. Stochastic resonance 



The deterministic dynamics of the species strongly depends on the value of the interaction parameter 0{t). We 

first analyze therefore the time evolution of /3{t) for different levels of the additive noise ap^'^. In the absence of noise 
((7/3 = 0) we have a periodical behaviour of (3{t) in the coexistence regime (see Fig. 2a). For low noise intensity with 
respect to the height of the potential barrier {a^ -C h) we obtain a quasi-deterministic trajectory of (i{t), slightly 
perturbed by the noise (Fig. 2b). When the noise intensity is comparable with the barrier height (cr^ ~ h), the typical 
picture of stochastic resonance phenomenon is obtained: the periodical behavior of the interaction parameter jumps 
between the two values f3 = 0.94 < 1 and /? = 1.04 > 1, which characterize the coexistence and the exclusion dynamical 
regimes (see Fig. (2c). A further increase of the noise intensity produces a loss of coherence and the dynamical behavior 
of P{t) is strongly controlled by the noise (Fig. 2d). 
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FIG. 2. Time evolution of the interaction parameter, with initial value /3(0) = 0.94, for different values of the additive noise 



intensity a p. (a) a-p = Q; (h) 
are 7 = 10"^ a;o/(27r) = 10"^. 



1.78 ■ 10""; (c) ag = 1.78 • lO"''; (d) ap = 1.78 • 10"^ The values of the other parameters 



The optimum coherent time behaviour of (i{t) (Fig. (2c) is obtained from the statistical synchronization relation 

Tk = To/2 (7) 
typical of the SR phenomenon. Here Tk is the Kramers time 

27r 

Tk = — = exp \2h/ag], (8) 

v/|C/"(0.99)|C/"(0.94) / '^J' ^ ' 

To = 2tt/ujo is the period of the driving force, [/"(0.99) and [/"(0.94) are the second derivative respectively calculated 
in the unstable and stable states of the potential. For Tg = 10^ we have ag = 1.78 • 10~^, the value of the noise 
intensity of Fig. 2c. 

Now we analyze the dynamics of the two species by fixing the additive noise intensity at the value cr^ = 1.78 • 10^"^, 
corresponding to a competition regime between the two species periodically switched from coexistence to exclusion, 
and by varying the intensity of the multiplicative colored noise. We obtain the temporal series of the two species for 
different values of both the multiplicative noise intensity a = ax = cry and the correlation time Tc- We do simulations 
in both regimes for the correlation time of the noise: (a) dt < Tc < To, and (b) Tc > Tg. In the weak correlated noise 
regime (a) no relevant modifications occur in the temporal series of the two species densities. The time evolution 
of the two species show an anticorrelated behavior with quasiperiodical oscillations (see Figs. 3c,4c) with a random 
inversion of the population that predominates over the other one, as in the white noise case In Figs. 3,4 we report 
these time behaviours at cr = 0, cr = 10~^^, a = 10~^, a = 10~^ for two different values of the correlation time 
Tc = 20, 2 • 10^. For < Tc <To = 10^ we obtain the same behaviour as in the presence of multiplicative white noise. 
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For Tc ~ To some modifications occur, particularly for Tc = 2 • 10"^ the temporal series of the two species densities 
show the anticorrelated behavior with quasiperiodical oscillations up to cr = 10~^ (see Fig. 4d), i. e. a delay in the 
coherent output of our ecosystem. This delay will manifest itself in the behaviour of the signal-to-noise ratio (SNR) 
as a function of the multiplicative noise intensity. 
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FIG. 3. Time evolution of both populations at different levels of the multiplicative noise for Tc = 20/ (a) a = 0; (b) a = 10~^^; 
(c) a = 10~*; (d) a = .The values of the parameters are /u = 1, a = 1, 7 = 10~^, a;o/(27r) = 10~^. The intensity of the 
additive noise is fixed at the value 0/3 — 1.78 ■ 10~^. The initial values are: for the two species a;(0) = j/(0) = 1, for the additive 
(white) noise /3(0) = 0.94, for the multiplicative (colored) noise Ci(0) = C2(0) = 0. 



In the strong correlated noise regime (b) a relevant delay of the coherent time behaviour of the two species is 
observed. The maximum of SNR is shifted towards higher values of the multiplicative noise intensity. This shift in a 
Log-Log scale grows faster than a linear function of the correlation time Tc. In Fig. 5,6 we report the temporal series 
for Tc = 2 • 10®, 2 • 10^ obtained for different values of the multiplicative noise intensity. 
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FIG. 4. Time evolution of both populations at different levels of the multiplicative noise for Tc = 2 - 10^; (a) a = 0; (b) 
a = 10~^^; (c) a — 10~* ; (d) a = 10~^. The values of the parameters, the intensity of the additive noise and the initial 
conditions are the same of Fig. 3. 
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FIG. 5. Time evolution of both populations at different levels of the multiplicative noise for Tc = 2 • 10 ; (a) a = 0; (b) 
a = 10~^^,- (c) a — 10^; (d) a = 10^. The values of the parameters, the intensity of the additive noise and the initial conditions 
are the same of Fig. 3. 



Comparing the temporal series of Figs. 3,4 and Figs. 5,6 we note that increasing the values of Tc the noise-induced 
effects are observed for higher levels of noise intensity. The coexistence regime and the correlated oscillations of both 
populations persist for a wider range of multiplicative noise intensities. The anticorrclatcd behavior with quasipcri- 
odical oscillations appear with very high noise intensity as the correlation time value of the multiplicative noise is 
strong enough. The lost of coherence in the time behaviours of the two species is observed at very high intensities of 
the multiplicative noise. Because of the high values of the multiplicative noise, one population extinguishes and the 
other one survives at a constant density after a transient dynamics (see Figs. 5d,6d). This dynamical behaviour is 
typical of an ecosystem in the presence of an absorbing barrier^. 

The periodicity of the noise-induced oscillations in the time behavior of the population densities shown in Figs. 3-6 
is the same of the driving periodic term of Eq.(6). 
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FIG. 6. Time evolution of both populations at different levels of the multiplicative noise for Tc = 2 • 10^; (a) a = 0; (b) 
a = 10^^^; (c) a — 10^; (d) a = 10^^. The values of the parameters, the intensity of the additive noise and the initial 

conditions are the same of Fig. 3. 
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FIG. 7. Log-Log plot of SNR as a function of noise intensity. The SNR is obtained for six different values of the correlation 

time: Tc = 2 ■ 10~^, Tc = 2 ■ lO""*^; Tc = 20; Tc = 2 ■ 10^, Tc = 2 ■ 10^; Tc — 2 ■ 10®. Moreover the signal-noise ratio for white 
Gaussian noise is reported. The SNR corresponds to the squared difference of population densities {x — y)^ . The values of the 
parameters and the initial condirtions are the same of those used to obtain the temporal series. 



This is the signature of stochastic resonance (SR) phenomenon. In order to manifest the presence of SR we analyze 

the squared difFcrcncc of population densities (x — y)^ for different values of Tc. In Fig. 7 the SNRs of this quantity are 
shown for Tc = 0, 2 • 10^^, 2 • 10^^, 20, 2 • lO'', 2 • 10®, 2 • 10^ as a function of the multiplicative noise intensity ct, by fixing 
the additive noise intensity at ap = 1.78 • 10~^. We calculate the SNR of {x — iff' because of the random periodical 
inversion of population densities. We note that dynamics of {x — y) is mainly affected by the multiplicative noise, as 
we can see from Eqs.(2). The SNRs have been obtained by performing 190 realizations of Eqs. (2). In each graph of 
Fig. 7 a maximum appears, whose position depends on the values of Tc, i.e. the most coherent response of the system 
is connected with both the intensity and the correlation time of the multiplicative noise. We see clearly the two 
dynamical regimes of weak correlated noise (a), with the first four values of Tc, and of strong correlated noise (b) with 
the last three values of t,,. In this regime the maximum of SNR is shifted towards higher values of multiplicative noise 
intensity as in previous theoretical and experimental studies®'^'^°. However some differences occur. Previous studies 
on the effect of colored noise on SR phenomenon showed that by increasing Tc the peak of the SNR shifts towards 
higher values of the noise amplitude and the maximum decreases with a broadening of the entire curve. The shift of 
the SR peak to larger noise intensities is due to the fact that colored noise suppresses exponentially the hopping rate 
with increasing noise color. In our model the colored noise is introduced in the multiplicative noise and not in the 
additive one as in usual bistable dynamical systems. The SR in the dynamics of interaction parameter /3 induces SR 
phenomenon in the dynamics of two competing populations^^. Our hopping rate in the first SR is not affected by the 
"color" of the multiplicative noise. However this noise produces the coherent response of the ecosystem and therefore 
produces the shift of the SNR peak. A more detailed study on the effect of colored noise on two noise sources used 
in our model will be the subject of future investigation. 



III. CONCLUSIONS 



We report a study on the role of the colored noise in the dynamics of two competing species. In our model we 
consider two Gaussian noise sources: a multiplicative colored noise and an additive white noise, which produces a 
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random interaction parameter between the species. The additive noise controls the switching between the coexistence 
and the exclusion dynamical regimes, the multiplicative noise is responsible for periodical oscillations of the two species 
densities, whose amplitude and coherence depend on the value of the correlation time Tc- When the interaction 
parameter is in SR regime a coherent time behavior of two species appears. This behavior is characterized by 
periodical oscillations and enhancement of the response of the system through the appearance of a stochastic resonance 
phenomenon in the dynamics of the two competing species. We then conclude that a stochastic resonance phenomenon 
(primary SR) for the interaction parameter /3 induces another effect of stochastic resonance (secondary SR) in the 
dynamics of the species in the presence of colored noise. We note that for Tc — > our results are consistent with 
that obtained for the case of white noise (sec Fig. 7). Moreover the coherent time behavior of our ecosystem and the 
SR phenomenon are shifted towards higher noise intensities, in agreement with previous theoretical and experimantal 
investigations^'^'^°. The presence of a colored noise assures that the densities of the two species at the time t depend 
on the values of the multiplicative noise at the previous times. This represents a more realistic condition and it can 
be useful to explain the time evolution of ecological species, whose dynamics is strongly affected by the environmental 
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